Improved qubit bifurcation readout in the straddling regime of circuit QED 
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We study bifurcation measurement of a multi-level superconducting qubit using a nonlinear res- 
onator biased in the straddling regime, where the resonator frequency sits between two qubit tran- 
sition frequencies. We find that high-fidelity bifurcation measurements are possible because of 
the enhanced qubit-state-dependent pull of the resonator frequency, the behavior of qubit-induced 
nonlinearities and the reduced Purcell decay rate of the qubit that can be realized in this regime. 
Numerical simulations find up to a threefold improvement in qubit readout fidelity when operat- 
ing in, rather than outside of, the straddling regime. High-fidelity measurements can be obtained 
at much smaller qubit-resonator couplings than current typical experimental realizations, reducing 
spectral crowding and potentially simplifying the implementation of multi-qubit devices. 

PACS numbers: 85.25.Cp, 74.78.Na, 03.67.Lx, 42.50.Lc, 42.65.Wi 



I. INTRODUCTION 

Circuit quantum electrodynamics (cQED), where su- 
perconducting qubits are coupled to transmission-line res- 
onators, constitute a promising architecture for the real- 
ization of a quantum information processor [l|, Q ■ Two 
criteria required for quantum computation are the imple- 
mentation, in a scalable way, of a universal set of gates 
and the ability to faithfully measure the qubit state Q. 
In this system, single qubit gates can be performed by 
sending microwave signals through the resonator close to 
the qubits' transition frequency, while two-qubit gates 
can be performed by tuning the qubits in and out of 
resonance. The increasing fidelity of one- Q and two- 
qubit [1-0] gates has allowed circuit QED to reach im- 
portant milestones, such as the implementation of two- 
and three-qubit quantum algorithms [8l-fl0| and the real- 
ization of more complex multi-qubit devices (Til ]. 

Qubit measurement in cQED is realized by driving 
the resonator close to its natural resonance frequency 
and by measuring the reflected or transmitted microwave 
signal. Recently, high-fidelity single-shot measurements 
have been achieved by using very large measurement 
drive powers [l2l - [l4| , by turning the resonator into a non- 
linear active device and using bifurcation to distinguish 
the qubit states I15l- 17j . or by using nearly quantum- 
limited amplifiers [18|. In these realizations, increasing 
the qubit-resonator coupling leads to larger variation of 
the resonator's parameters with the qubit state, resulting 
in high measurement fidelity. In the same way, increasing 
this coupling also typically reduces the gate-time of two- 
qubit operations. However, stronger coupling can also 
reduce the on/off ratio of logical gates, causes spectral 
crowding and reduces the qubit lifetime through spon- 
taneous emission via the resonator, also known as the 
Purcell effect. 

In this paper, we take a different approach and show 
that it is possible to implement high-fidelity single-shot 
measurements of a superconducting qubit using relatively 
small qubit-resonator coupling strengths — of the or- 



der of 10 MHz — than in many recent experiments. 
To achieve this, we use the weakly anharmonic multi- 
level structure relevant for most superconducting qubits 
and take advantage of the so-called straddling regime 
where the resonator frequency sits between two qubit 
transitions [l9T ] . This regime shows enhanced qubit-state- 
dependent pull of the resonator frequency, enhanced 
qubit-induced resonators and reduced Purcell decay rate. 
We show that these three characteristics combine to im- 
prove bifurcation measurements of the qubit state. In nu- 
merical simulations of qubit readout, we find error proba- 
bilities three times smaller inside with respect to outside 
of the straddling regime. Even without thorough explo- 
ration of the available parameter space, we find measure- 
ment fidelities of 98%. 

The paper is organized as follows. In section [Hi we 
first introduce the Hamiltonians modeling a nonlinear res- 
onator, required for bifurcating measurements, coupled 
to a multi-level qubit. Then, in section IIII1 we review 
the principle of bifurcation measurements and highlight 
the important differences between two-level and multi- 
level qubits in this respect. In section IIV1 we derive an 
effective dispersive Hamiltonian valid in the straddling 
regime. Finally, we compare in section fVl the parameters 
calculated with our model to parameters extracted from 
exact diagonalization of the qubit-resonator Hamiltonian. 
We then examine the specifics of bifurcation in the strad- 
dling regime, extract measurement fidelities from numer- 
ical simulations and discuss other advantages of working 
in this regime. 



II. MODEL 

As mentioned above, many superconducting qubits 
have a relatively small anharmonicity and are therefore 
described by M-level systems with M > 2 rather than 
by two-level systems |19h22| |. We consider such a qubit 
coupled to a Kerr nonlinear resonator, which could be 
realized for example by an LC-circuit with a Joseph- 
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son junction II 5[ or a stripline resonator with one [17[ 
or many [23l . |24| | embedded Josephson junctions making 
it nonlinear. The qubit-resonator system can be mod- 
eled with the many-level version of the Jaynes-Cummings 
Hamiltonian 

H s = H q + H r + Hi, (2.1) 

where (h = 1) 

M-l 

H q =Y. ^i,i=^, (2.2) 

i=0 

is the qubit Hamiltonian, 

K 

H r = bj r a) a + — aWaa, (2-3) 
is the nonlinear resonator Hamiltonian [25j . and 

M-2 

Hl =H ft(o t n j , i+ i + an i+1 ,i), (2.4) 

is the interaction Hamiltonian and where Iljj = |i) (j| 
with the qubit eigenstates. In these expressions, uji 
is the frequency associated to the qubit eigenstatc uj r 
is the bare resonator frequency (at low powers), K is the 
Kerr constant, and gi the qubit-resonator coupling con- 
stants. We have also introduced the short-hand notation 

M-l 

n, = •''•"<•• ( 2 - 5 ) 

i=0 

where a; is a scalar taking different values Xi associated 
to the different qubit states \i) . This notation is used 
throughout this paper. Finally, in the qubit-resonator in- 
teraction term, we have made the standard rotating- wave 
approximation (RWA) and also assumed that transition 
between states \i) <-> \j) are suppressed for \i— j\ ^ 1 [l9j |. 

Measurement of the qubit is realized by driving the 
resonator with a tone of amplitude and frequency uid- 
This is modeled by the drive Hamiltonian 

fl* = ed^'V + e^'a), (2.6) 

leading to the total Hamiltonian 

H = H S + H d . (2.7) 

III. BASICS OF BIFURCATION 
MEASUREMENTS 

The description of the Kerr nonlinear resonator (KNR) 
is simplified by introducing the reduced detuning fre- 
quency VL = 2(co r — LUd)/n J26|. As illustrated in Fig. [TJ 
the steady-state response of the KNR can vary drasti- 
cally whether the reduced detuning D, is larger or smaller 
(in absolute value) than a critical detuning ile = V3- 
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FIG. 1. (Color online) (a) Response (amplitude of the field) 
as a function of the reduced detuning and for increasing drive 
amplitude e p . The back-bending of the response reflects the 
choice K < 0, as is usually the case in circuit QED (b,c) Sta- 
bility diagram of a Kerr nonlinear resonator in the absence (b) 
or presence (c) of a qubit. Inset : Time-dependent enveloppe 
of a sample-and-hold readout. 

For |f2/f2c| < 1, the resonator response is single valued, 
with as shown in Fig. [T] (a), a response that is stiffened 
compared to the usual Lorentzian line shape. Close, but 
below, the critical point, the resonator can then be used 
as a parametric amplifier for small signals (23|. On the 
other hand, for |f2/f2c7| > 1 the resonator is in the so- 
called bifurcation amplification regime (BA) where it is 
bistable for a range of drive amplitudes e p . If e p is ramped 
up starting from zero, the resonator's response will bifur- 
cate from a low (L) to a high (H) oscillation amplitude 
dynamical state at a critical amplitude en- If the drive 
amplitude is then reduced, the resonator stays in the 
state H until the drive amplitude becomes lower than 
a second threshold ei- The associated stability diagram 
is illustrated in Fig. Q] (b). 

As was already experimentally demonstrated, in the 
BA regime, the KNR can be used as a sample-and-hold 
detector of a qubit [H, [13, [H HI] . Indeed, as for most 
quantum information related tasks, qubit readout is re- 
alized in the dispersive regime where \gi\ <§; \u> r — w-s+i^l. 
In this situation, the system Hamiltonian H s is well ap- 
proximated by the effective Hamiltonian [29| 

K 

Hd pa {uj r + Hs)a^a + —aJa'aa + n w . (3-1) 
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As can be seen from the coefficient of a^a, in this regime, 
the presence of the qubit results in a shift of the res- 
onator frequency by a qubit-state dependent quantity Si. 
This dispersive cavity pull, whose value Si will be dis- 
cussed below, results in different thresholds e^j and en,i 
depending on the qubit states. This is schematized for 
the first two qubit states {|0) , |1)} by the red and the 
blue lines in Fig. [T] (c). 

Starting from zero, increasing the drive amplitude ed 
until en.o < £d < e H,i will then result in a high amplitude 
of the cavity field if the qubit is in its ground state and 
a low amplitude if it is in its excited state. This range 
is represented by the gray shaded area in Fig. [1] (c). If 
the drive amplitude is then reduced below en,o, but stays 
above cl,i/o [ see inset of Fig. Q] (c)] , both resulting states 
are stable and the qubit state has been mapped into the 
dynamical state of the resonator. Since these dynamical 
state are stable, it is possible to accumulate the output 
signal for a time longer than the qubit relaxation time 
T\. The measurement fidelity can then be optimized by 
varying the sampling time t s , the height of the plateau 
5e m and the steepness of the ramp up a 

nana nam. 

In practice, the readout fidelity is limited b y q ubit re- 
laxation during or before the sample phase [17j . when 
the resonator has not bifurcated yet. The speed at which 
the sampling can be made is limited by the resonator's 
decay rate k. Indeed, ramping up the drive much faster 
than 1/k will produce large ringing oscillations in the 
field amplitude and which can result in false positives 
or negatives. This results in a reduced measurement fi- 
delity. Increasing k therefore implies smaller transients 
and hence faster measurement. However, increasing k 
too much can also yield a lower measurement fidelity. In- 
deed, in the limit where k is much larger than the differ- 
ence between the qubit-state dependent resonator pulls 
X = Si — So, both qubit states are indisthinguishable. 
Moreover, increasing k also increases t he q ubit 's Purcell 
decay rate 7 K ~ ngf /(uii+i — uii — w r ) 2 |3Cll. [3l| which ul- 
timately limits the qubit relaxation time T\ . Ideally, one 
would like to increase both x an d without increasing 
the Purcell rate. 



A. Two-level systems 

The qubit-state dependent resonator shift Si discussed 
above depends on the coupling <?i and the qubit-resonator 
detuning Aj jr = — u>i — oj r . For a two-level qubit, it 
takes the simple form [l[ 



Si 



2LS 

/0 



± 



A() ,T 



(3.2) 



corresponding to symmetric displacement of the cavity 
frequency around its bare frequency u r . The difference 
between the pulled resonator frequency for the qubit 
states |0) and |1) is therefore x — &i ~ So — %9o/^o,r- 
This cavity pull can be of the order of a few tens of MHz, 
while staying in the dispersive regime, with the typical 



(a) 

Co>2 
UJq 



Qubit 
frequencies 









Wio 



(b) 



Outside 
Straddling 



(c) 



Inside 
Straddling 



U21- 



FIG. 2. (Color online) (a) Absolute and transition frequencies 
of the first three eigenstates of a qubit with negative anhar- 
monicity such as the transmon. (b) Example of frequency 
diagram where the resonator frequency (dashed lines) is out- 
side of a straddling regime, (c) Example of energy diagram 
where the resonator frequency (full black line) is inside of a 
straddling regime. 



values go/2-K ~ 100 - 200 MHz and A , r /27r - 1 - 2 GHz. 
Such couplings have been achieved with transmons and 
flux qubits [32|, [33|. For a two-level qubit, increasing the 
coupling <?o increases x, but also increases j K by the same 
amount. For Purcell-limited qubits, this negates the gain 
of this strategy. 



B. Multi-level systems 

For multi-level systems, the shifts are changed by the 
presence of additional levels, and the symetry around the 
bare resonator frequency is broken. Indeed, the frequency 
shift is given by [19( 



:MLS 



2 

9i-i 



(3-3) 



As illustrated in Fig.f2](b), in most experiments 0,0, HH 
the qubits are biased such that the resonator frequency 
sits above, or below, all of the qubit transition frequen- 
cies. This results in a pull of the resonator frequency 
X = 2<?o/Ao,r — Si/Ai.n reduced compared to that 
of a purely two-level system. In the limit where the 
multi-level system tends toward a harmonic oscillator, 
Ai jf . — > Ao, r and <?i — > \/2.9o such that this pull vanishes. 
The reduction in the pull can be compensated with larger 
couplings gi achieved for example with transmons [20| . 
However, as stated above, increasing the qubit-resonator 
coupling also increases the resonator-mediated Purcell de- 
cay 31 1 and dressed dephasing (33 - [3r| . This dependence 
on k and gi of both resonator-mediated qubit decay and 
measurement speed ultimately limits the achievable mea- 
surement fidelity. 

One way to increase the dispersive shifts x without 
increasing the coupling is to work in the so-called strad- 
dling regime [l9[. In this regime, illustrated in Fig.[2](c), 
the detunings Ai !r and Ao ir are of opposite signs. As a 
result, instead of canceling each other, the two terms in 
Eq. (|3.3[) add up, yielding a significantly enhanced value 
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of x- Since this improvement is obtained without increas- 
ing <7j, it does not increase the Purcell rate. Moreover, 
as we will show in the next section, this regime also in- 
creases qubit- induced nonlinearities [l3j . something that 
we will exploit below to improve bifurcation readouts. 



IV. DISPERSIVE MODEL IN THE 
STRADDLING REGIME 

Following the approach of Ref. [2!| , we use a polaron 
transformation [37H39{ followed by a dispersive transfor- 
mation [U, |4(j to approximately diagonalize the Hamil- 
tonian of Eq. (|2.7[) . Doing the transformations in this or- 
der (polaron followed by dispersive), allows to correctly 
model the ac-Stark shift caused by a drive detuned from 
the resonator frequency (29j . However, since we are inter- 
rested in the straddling regime, one more transformation 
must to be done in order to diagonalize an effective two- 
photon process that is important only in the straddling 
regime. This is done in Appendix [A"l and yields the effec- 
tive diagonal Hamiltonian 

H' = IL U , + [u' r (a) + n s(Q) ] a t a , (4.1) 

where u)' r (a) is the Kerr-shifted resonator frequency 

uj' r (a) = uj r + 2K\a\ 2 , (4.2) 

with a = (n a ) the resonator mean field and 

to'i =u>i + S.lal 2 + K t \a\ 4 + L^a), (4.3) 

are the renormalized effective qubit frequencies. There, 
we have defined 



-il 2 ) 



— (Xi — Xi_i), 
-SidAil 2 + |Ai 
-±(3Xi +1 |Ai| 2 -JQ|A i+1 | 2 ) 



(4.4a) 
(4.4b) 



X 



(2) 



X 



(2) 



the linear and quadratic ac-Stark shift coefficients with 
A* = —gi/&i,d, Xi = -gih-i, and where A^ d = lu 1+1 - 
LUi — uod is the detuning between the qubit transition i 
and the drive d. The last line of K, comes from the diag- 
onalization of an effective two-photon transition process 
that is large only in the straddling regime. This con- 



tributes the last two terms of K,- with X 



A 



(2) 



-- -g\ 7(A; +M + A M ) and where 

g {2) =AiA i+1 (A i+lid -A iid ). 
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(2) A (2) 



(4.5) 



We note that, when compared with the results of Ref. [13J, 
the detunings Aj are defined with respect to the drive 
frequency, and not the resonator frequency. In addition, 
in Ref. 1 1 31 ] , the dispersive transformation was done with 
respect to the field operator a rather than to the classical 



field a. Because of this choice, the quadratic term Sj in 
Ref. [IH contains correction which accounts for a specific 
choice of ordering for the ladder operators in the quartic 
term. Here, since it is the classical field that is considered, 
there are no such corrections (i.e. aa* = a* a). 

Finally, in Eq. (|4.3[) we have also defined the Lamb 
shift 



L» (a) = 



.97 



where uj'/ is given by 



^lal 4 . 



(4.6) 



(4.7) 



Using this definition, the cavity pull II,g( Q ) in the effective 
Hamiltonian Eq. (|4.1j) can be expressed in a compact way 
using Si(a) = — [Lj+i(a) — Li(a)}. We note that while 
the ac-Stark shift coefficients Si and Kj depend on the 
qubit-drive detuning, the Lamb-shift Li depends on the 
detuning between the ac-Stark shifted qubit and Kerr- 
shifted resonator (29[. Finally, the steady-state qubit- 
state-dependent cavity field a.i is given by the solution 
of 



-e d = (u r -Wd-i%)a i +K\a i \ 2 a i + S, + 



4 



(4.8) 



In Fig. |3l we compare the above analytical expressions 
for Sj and Kj to numerical results. These quantities 
are found numerically by fitting a quadratic polynomial 
to the resonator frequency for the qubit state \i) and 
in the presence of n photons, uj r ^(n) = JSj jtl +i — E^ n . 
The energy E^ n is found numerically by diagonalizing 
the undriven qubit-resonator Hamiltonian H s and tak- 
ing K = 0. We then associate Ei in to the energy of the 
cigenstate closest to the bare qubit-resonator state \i,n). 
The parameters, given in the caption of Fig. [3j arc typical 
to transmon qubits [l9| , but with a smaller than typical 
coupling gio/2ir = 13.5 MHz. We show the analytical [(a) 
and (c)] and numerical [(b) and (d)] values of Si [(a) and 
(b)] and K, [(c) and (d)] for the ground state i = (full 
black lines) and first qubit excited states i = 1 (dashed 
green lines). We find quantitative agreement, except at 
the qubit-resonator resonances and at the two-photon res- 
onances (identified by divergences). We finally show in 
Fig. [3] (e) the Purcell decay rate 7 K i of level i assuming 
uj r = cJd and k/2tt = 5 MHz. 

Two operating points, designated by A and B and iden- 
tified by the vertical full black lines and dashed grey lines 
respectively, are illustrated on Fig. [3] These particular 
points have been chosen because, while A lies in the strad- 
dling regime and B is outside of that regime, the cavity 
pull |x| = |§i — So| is identical in both cases. In the 
next section, we will show that working in the straddling 
regime is advantageous for qubit readout. Since the cav- 
ity pull is the same at both A and B, improvement in the 
measurement will be due to qubit-induced nonlinearities 
Kj or variation in the Purcell decay rate. 
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FIG. 3. (Color online) [(a) and (b)] Linear and [(c) and 
(d)] quadratic ac-Stark shifts for the ground (full black 
lines) and the excited (green dashed lines) states of a trans- 
mon qubit with charging energy Ec = 300 MHz, 

Josephson energy Ej = 25 GHz, coupling at zero flux 
gio/2-K = 15 MHz tuned such that uj w /2n = 6 GHz. This 
yields (u , ui, w 2 , w 3 )/2tt « (0,6,11.7,16.9,21.8) GHz and 
(5io,P2i,532,543)/2tt = (13.5,18.5,21.8,24.1) MHz at the op- 
erating point. Panels (a) and (c) correspond to the analytical 
equations Eq. (|4,4|1 . while panels (b) and (d) are extracted 
numerically as described in the text. Panel (e) shows the Pur- 
cell decay rate 7«i = Kgt/A? r assuming k/2-r = 5 MHz and 
u r — ujd for i = (full black line), i = 1 (dashed green line) 
and i = 2 (dotted blue line) . The black line, corresponding to 
the |1) — > |0) transition, is the relevant one for qubit operation. 
The shaded area corresponds to the straddling regime. The 
full vertical black line and the dashed vertical grey line cor- 
respond to the two operating points A (ud/2n — 5720 MHz) 
and B (uj d /2-K = 6044 MHz) discussed in Sees. HVl and M 



The qubit-induced nonlinearities Kj are plotted in pan- 
els (c) and (d). Comparing panels (a) and (c), we note 
a major difference between the operating points A and 
B. At B, the sign of Kj is opposite to that of §i for both 
? = (full black lines) and i = 1 (dashed green lines). 
This sign difference corresponds to a cavity pull that is 
decreasing when the number of photons increases. On 
the other hand, at point A, the sign of Kq is the same 
as that of § ■ Therefore, we expect that the cavity pull 
at point A will not decrease as much as at point B with 
increasing photon number [i~3j . Moreover, we can see in 
panel (e) that the Purcell rate for the transition |1) — ^ |0) 



(full black line) is much larger at point B than at point 
A. 

One would expect that these two effects — a cavity pull 
that reduces less with increasing number of photon and 
a reduced Purcell decay rate — lead to better qubit mea- 
surement at operating point A than B. In the next sec- 
tion, we show numerically that this expectation holds for 
a Kerr resonator operated close to its bifurcation point. 
This is done by first calculating the steady-state photon 
number associated to both qubit states. We then sim- 
ulate the complete dynamics corresponding to a qubit 
under measurement with the microwave pulse typically 
used in bifurcating readouts [13, H3, E| and which is de- 
signed to make the resonator latch in its H state for one 
of the qubit state. From these simulations, we extract 
the expected measurement fidelity and show that better 
results are indeed obtained at operating point A than B. 



V. IMPROVING BIFURCATION 
MEASUREMENTS IN THE STRADDLING 
REGIME 

Bifurcation measurements rely on the critical drive am- 
plitude en,i — at which the resonator bifurcate to its high 
state H — being different for each qubit state i. As illus- 
trated in the inset of Fig. [1] (c), in bifurcation measure- 
ments the measurement drive amplitude is increased to 
a value between these two critical amplitudes. However, 
the bifurcation process being probabilistic, the resonator 
can still bifurcates from the L to the H state even if the 
drive amplitude is (slightly) lower than en- This yields 
errors in the measurement and a reduced measurement 
fidelity. We therefore expect the measurement fidelity to 
increase with Ae# = \en.o ~ e H,i\ and so, in other words, 
with cavity pull. In addition, one expects that a larger 
separation of the thresholds protects the measurement 
better against ringing in the resonator's response which, 
close to £h , may lead to unwanted bifurcation. For these 
reasons, we expect that the operating point A, at which 
the cavity pull should remain larger on a wider range of 
measurement power, to be better for measurement than 
point B. 

Below, we first calculate the steady-state response of 
the resonator in section IV Al We then compute the mea- 
surement fidelity for a pulsed measurement in section lVBI 
Finally, we discuss other advantages of working in the 
straddling regime in section IV CI 



A. Steady-state response 

We simulate the evolution of the state p starting with 
the resonator in the vacuum and with the qubit either 
in the eigenstate |0) or |1). We first focus on a drive of 
constant amplitude e^, without intrinsic qubit relaxation 
or dephasing. By looking at the resonator's steady-state 
response, with this simulation, we want to show that the 
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FIG. 4. (Color online) Numerically computed average num- 
ber of photons rii for a qubit initialized in state \i) with 
i = (a) and i = 1 (b) without qubit relaxation or de- 
phasing. The evolution is computed according to Eq. (|5.ip . 
Qubit parameters are given in the caption of Fig. [3] Res- 
onator parameters are (n,K)/2ir = (5,-0.4) MHz and the 
resonator frequency is adjusted to keep the resonator- drive 
detuning (lj t — (jjd)/2ir = 15 MHz such that the reduced de- 
tuning 0,/Q.c ~ 3.5, well in the bifurcation regime. Dashed 
lines represent the two operating points A and B (see caption 
of Fig. [3|. Dots indicates the bifurcation thresholds en,i (red 
for i = 0, blue for i — 1). 



distance between the bifurcation thresholds en,o and 
is indeed larger at operating point A than B. The evolu- 
tion is governed by the master equation 



P 



-i [H,p] + nT>[a]p, 



(5.1) 



-a^ap- 



with the Linblad-form dissipator T>[a] = i(2apa^- 
pala). After a time long compared with 1/k, we compute 
the average number of photon for the qubit initially 
in state i £ {0, 1}. 

This quantity is plotted in Fig. @] as a function of the 
drive frequency Ud and amplitude ed for the qubit initially 
in its ground (a) or excited (b) state. In both cases, two 
regions corresponding respectively to the resonator being 
in the L state (dark blue, n., < 10 photons) or in the H 
state (light red, ~ 50 photons) can be identified. The 
border between these two regions (white) corresponds to 
the critical drive amplitude €H,i, at which the photon 
population sharply goes from rii ~ 15 to rii ~ 50. When 
comparing these results to the dispersive shifts illustrated 
in Fig. |31 we can see that sharp changes in Si and Kj trans- 
late into sharp changes in the bifurcation amplitudes ch,i- 
For example, both So and Si [Fig. [3] (a)] change sign at 
uid/'Z'K = 6 GHz, which translate in a sharp change in 
both en,i around that frequency. Moreover, §i changes 
sign at Wd/27r ps 5.7 GHz, while Si does not. As a result, 
as can be seen in Fig- HI only ch,i changes significantly at 
that frequency. Finally, variations in Kj are also visible, 
for example as the feature in eh,o at Ud/2-K ps 5.85 GHz 
corresponding to the change of sign in Ko at that same 
frequency. 

The operating points A and B are illustrated in Fig. @] 
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FIG. 5. (Color online) Typical Q function Q(a) of a resonator 
when driven close to its bifurcation threshold en- 



by the horizontal full black lines and dashed gray lines 
respectively. The thresholds es,i at these two points are 
identified by full circles (red for €h,o and blue for £h,i)- 
As expected from the above arguments, the separation 
Aen = \ch,o — e H,i\ is larger at A than at B. For the cho- 
sen parameters, we find Aejj /2tt ~ 10 MHz at A while we 
find Aen/^n ~ 5 MHz at B. We note that €h,i > (H.a at 
point A while ch,i < £h,o at point B. This simply changes 
which resonator state — of L and H — is associated with 
each qubit state. 



B. Pulsed measurement fidelity 

In order to quantify by how much an actual measure- 
ment can be improved by working at operating point A 
- inside the straddling regime — rather than at B 
outside of the straddling regime — we numerically simu- 
lated a bifurcation measurement with a sample-and-hold 
shaped pulse as illustrated in the inset of Fig.[T](c). We 
recall that, to our knowledge, all experiments with bi- 
furcation measurements have been made outside of the 
straddling regime so far. 

To be more realistic, we performed numerical inte- 
gration of master equation Eq. (|5.ip including qubit 
dissipation modeled using the Lindblad-form term 



7E,=o V 
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j-nj^_|_i p. Here, 7 is the decay rate of the 

first qubit transition and the factor gi/go is included to 
take into account the variation of the qubit decay rate 
with increasing i [29| . Pure dephasing is not included 
since recent devices tend to have very low pure dephas- 
ing rates [13, [H|- Including this effect would possibly 
affect the QND character of the readout due to dressed 
dephasing [3J-|3g , but the extent of this effect has yet to 
be measured experimentally. 

At the end of the hold time, the Q function of the 
resonator Q{a) = (a \p\ a) /ir is computed. A typical Q 
function near the bifurcation threshold en is represented 
in Fig. [5] It shows two well-separated smooth peaks cor- 
responding to the L and H states of the resonator. The 
switching probability is extracted from the weight of the 
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FIG. 6. (Color online) Error probability for the outcome of a 
bifurcation measurement versus sampling time t a . The mea- 
surement pulsed is illustrated in the inset of Fig. Q](c). Full 
black lines correspond to operating point A, in the straddling 
regime and dashed grey lines correspond to operating point 
B, outside of the straddling regime. The parameters are the 
same as in Fig. [4] The simulations are realized for a qubit 
with intrinsic relaxation times (a) Ti = oo, (b) 800 ns and (c) 
300 ns. In (b) and (c), the dashed vertical red lines indicate 
the minimum of the curves for the two operating points. Full 
green bars indicate the gain in measurement time (horizontal) 
or measurement fidelity (vertical) between the two operating 
points. 



peak that is the farthest away from the origin. From the 
switching probabilities, the worst-case error probability 



Perror = max P(j)\l), 

{i,j'}e{o,i} 



(5.2) 



can be computed and where P{j\i) is the probability of 
assigning the measurement to the qubit state given 
that the qubit was initially in \i). This numerical pro- 
cedure was previously tested against experimental single- 
shot bifurcation measurement of a transmon qubit [17[ 
and found an identical measurement fidelity, within a 
margin of 2% 0. 

We show in Fig. [5] the worst-case error probability as a 
function of the sampling time t s for three different qubit 
relaxation times T\ . These results have been obtained by 
minimizing the error probability with respect to a and 
Sed [see inset of Fig Q] (c) for definitions] . Comparing 
panels (a), (b) and (c), we see that P crloI increases as 
the qubit relaxation time decreases, which is expected 
because of the increased odds of the qubit relaxing before 
the resonator switches from L to H . 

We now compare the results inside (full black lines, 
operating point A) and outside (dashed grey lines, oper- 
ating point B) of the straddling regime. We first observe 
that for short sampling times t s , the error probability 
is always lower for operating point A than B. Since the 
low-photon cavity pulls \ were chosen to be the same for 
both points, this improvement is due both to the sign 
and amplitude of the Kerr terms K, and to the reduced 
Purcell decay as explained in section IIV1 The situation 
is however reversed for larger t s where point B is supe- 



rior. As illustrated in Fig. [4] this is because the resonator 
switches at a lower power for the ground state than for 
the excited state at point A, while the opposite is true 
for point B. This implies that qubit relaxation induces 
resonator switching (ie false positives) at point A, but 
not at point B. We note that the situation would be re- 
versed for a qubit with a positive anharmonicity such as 
the low- impedance flux qubit 22j, increasing further the 
advantage of working in the straddling regime. 

Overall, we find that operating within the straddling 
regime always allows to reach lower error probabilities 
with a sampling time t s always as short, or shorter, than 
outside of the straddling regime. When operating in the 
straddling regime, the error probability is up to 3 times 
smaller than outside. Finally, the absolute improvement 
is better for qubits with shorter lifetimes, but as expected 
the best fidelity is found for qubits with longer lifetimes. 



C. Other advantages 

The above improvement in readout fidelity has been ob- 
tained by working with a qubit-resonator coupling g that 
is more than an order of magnitude smaller than current 
experimental realizations. Lower coupling however leads 
to slower two-qubit gates when these rely on qubit-qubit 
interactions mediated by the same resonator mode that 
is used for readout. This problem can be sidestepped 
by either taking advanta ge o f different modes for read- 
out and two-qubit gates [43J or, as recently experimen- 
tally realized, using direct capacitive coupling between 
the qubits Hi HI]. 

With the above problem avoided, working with weaker 
coupling g can be advantageous in other ways than the 
more efficient readout studied here. For example, it al- 
lows to greatly reduce Purcell decay by biasing the qubit 
away from a resonator resonance when it is not being 
measured. With a reduction by a factor of 10 of the 
coupling, a reduction by a factor of 100 of the Purcell 
decay rate can be obtained for the same detuning and 
cavity damping [see full black line in Fig. (He)]. At the 
time of measurement, the qubit-resonator detuning can 
be adjusted such as to reach the straddling regime. This 
can be done by changing the flux in th e q ubit loop or 
by using a tunable resonator (or both) [461 ] . Moving in 
and out of the straddling regime in this way necessarily 
means going through a qubit-resonator crossing. With a 
large coupling g, the associated (and unwanted) Landau- 
Zener-Stueckelberg transitions can be correspondingly 
large [47]] • This probably is however greatly reduced 
when working with small couplings. Indeed, assum- 
ing a frequency-tuning speed of v = 2tt ■ 1 GHz/1 ns, 
one finds the probability of unwanted transition P = 
1 _ e -^9?>/v = 0.7% for a coupling g/2iT = 13.5 MHz, 
while the same probability is ~ 10% for g/2ir = 50 MHz 
and - 30% for g/2n = 100 MHz. 

Smaller coupling strengths can also help in reducing 
spectral crowding in the presence of multiple qubits cou- 
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pled to a single resonator. Indeed, even if the qubit- 
qubit interaction mediated by virtual excitations of the 
resonator is not actively used for logical gates, it is al- 
ways present and can lead to errors. The rate of this 
interaction can be reduced by increasing the qubit-qubit 
detuning by an amount that is large with respect to the 
coupling g. With large g and multiple qubit, the avail- 
able spectral range (typically from ~ 4 to 15 GHz) is 
rapidly occupied and only a few qubits can be coupled 
to the same resonator without having to deal with un- 
wanted two-qubit gates. Using the straddling regime to 
increase the measurement fidelity with smaller coupling 
addresses this problem and does not require advanced 
circuit des: 



tins pro 
igns [IH 



VI. CONCLUSION 

We have studied the measurement of a multi-level su- 
perconducting qubit using bifurcation of a Kerr nonlinear 
resonator and by exploiting the straddling regime. The 
method is applicable to any qubit with a weakly anhar- 
monic multi-level structure with only nearest-level tran- 
sitions, but could be generalized to more complex struc- 
tures and couplings. As we have shown, working in the 
straddling regimes allows larger qubit-state-dependent 
pulls of the resonator frequency for a given coupling or, 
equivalently, the same pull for smaller couplings. While 
outside of the straddling regime, the resonator frequency 
shift is reduced at higher photon numbers [34| . we show 
that, inside the straddling regime, it is possible to find 
operating points where this reduction is minimized. We 
also show that the Purcell decay rate can be much smaller 
for a given cavity pull inside the straddling regime. Com- 
bined, these two effects lead to an increased fidelity for bi- 
furcation measurements and we find an error probability 
up to three times smaller inside than outside of the strad- 
dling regime for a sampling time that can be more than 
250 ns shorter. We find measurement fidelities 1 — P e rror 
larger than 98% with a qubit-resonator coupling as small 
as 13.5 MHz with realistic system parameters. 

The method presented in this paper has also the ad- 
vantage of reducing spectral crowding in multiple-qubit 
systems. It does that without requiring complex circuits 
and allows to effectively remove Purcell decay when the 
qubits are not being measured. 
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Appendix A: Dispersive transformation of the 
two-photon terms 

In this Appendix, we follow Ref. [1^] to diagonalizc 
the Hamiltonian (|2.7|) as well as a two-photon transition 
term that can be large only in the straddling regime. To 
do so, we first apply a polaron transformation [37H39j[ 



M-l 



P = ]T Hi,iD(ai), 



(Al) 



where D(a) is a displacement transformation [48| 

D(a) = exp [aa^ — a*a] , (A2) 

that displaces the resonator field operator a — > a + a,. 
The result of the polaron transformation on a is therefore 
a — > a + n a , where H a is defined according to Eq. (|2.5[) . 
We follow this polaron transformation by a dispersive 
transformation of the classical detuned drive on the qubit 



exp 



A/-2 



J2 Cn,,, + i - &n 



i+l, 



(A3) 



where £j is a classical analo gue of the operator A,a^ in the 
dispersive transformation [34l . [40j ] . Applying these two 
transformations on the Hamiltonian (|2.7[) and choosing 
on according to Eq. (14. 8|) and 



yields the Hamiltonian [2£ 

M-i 



-giCti 



->i+l,i 



(A4) 



H 



E 

Af-3 

E 



H 



r + uj' r (a)a' a 
^'^ILi+a + h.c. 



(A5) 



where the dispersive transformation has been performed 
to fourth order and g\ , given at Eq. (|4. 5[) , is an effective 
coupling due to two-photon transitions. In the above 
Hamiltonian, we have defined the ac-Stark shifted qubit 
frequencies 



uj'I = LJi + : 



d 1} |al 4 , 



(A6) 



where a = (IL a ), a, is given by the solution of Eq. (|4.8[) . 
S{ is given at Eq. (|4.4a|) . while is given by the first 
three lines of Eq. (|4.4b[) , and the Kerr-shifted resonator 
frequency 



ui' r (a) 



2K\a\ 2 . 



(A7) 



We note that the second line of H" is not diago- 
nal. In Ref. [2{|, this term was dropped assuming that 

(2) 

g\ was small and that \2ujd — w; + 2 + was large 
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enough to do a RWA. Here however, since we are in- 
terested in the straddling regime, the same can not be 
done. Indeed, if, for example, the drive frequency is 
u c i = (k>2 — 0Ji)/2 + (oji — Wo)/2, which falls directly 
in the middle of a straddling regime, the second line of 
H" is resonant and a two-photon transition from |0) to 
|1) is driven. Moreover, since Aj+i^ and A^ have the 

(2) 

same sign, the coupling g\ can be large. We can how- 
ever approximately diagonalize this term using a second 
transformation of the form 



D (2) 



exp 



M-3 



(2) 



n 



i,i+2 



t(2) n 



(A8) 



Applying this transformation on Eq. (|A_5|) and choosing 



t (2) 



(2) 



A i+1 , d + A 



i , d 



(A9) 



yields a correction to the Kerr shift, giving Eq. (|4.4bl) . 
Applying a final dispersive transformation on H" in or- 
der to diagonalize the quantum interaction Hi yields the 
diagonalized Hamiltonian (|4.1I) . 
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